Non-adiabatic electron pumping: maximal current with minimal noise 
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The noise properties of pump currents through an open double quantum dot setup with non- 
adiabatic ac driving are investigated. Driving frequencies close to the internal resonances of the 
double dot-system mark the optimal working points at which the pump current assumes a maxi- 
mum while its noise power possesses a remarkably low minimum. A rotating-wave approximation 
provides analytical expressions for the current and its noise power and allows to optimize the noise 
characteristics. The analytical results are compared to numerical results from a Floquet transport 
theory. 
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In mesoscopic conductors, a cyclic adiabatic change 
of the parameters can induce a pump current, i.e., a 
non-vanishing dc current flowing even in the absence of 
any external bias voltage [1~3]. For adiabatic quantum 
pumps [4-8], the transfered charge per cycle is deter- 
mined by the area enclosed in parameter space during the 
cyclic evolution [4, 5]. This implies that the resulting cur- 
rent is proportional to the driving frequency and, thus, 
suggests that non-adiabatic electron pumping is more ef- 
fective. For practical applications, it is also desirable to 
operate the quantum pump in a low-noise regime. It 
has been found that adiabatic pumps can be practically 
noiseless [9]. This happens, however, on the expense of 
acquiring a small or even vanishing current [10]. There- 
fore, the question arises whether it is possible to boost the 
pump current by increasing the driving frequency while 
keeping the noise level very low. 

Non-adiabatic electron pumping can be achieved ex- 
perimentally with double quantum dots under the in- 
fluence of microwave radiation [11-14]. In this letter 
we study the transport properties in this non-adiabatic 
regime. Our main aim is to find ideal parameter regimes 
in which a large pump current is associated with low 
current noise. For the optimization of the system pa- 
rameters, it is beneficial to obtain, besides a numerical 
solution, also analytic expressions for the transport quan- 
tifiers. Therefore, within a rotating-wave approximation 
(RWA) , we map the driven transport problem to a static 
one which we solve analytically. In doing so, a particular 
challenge represents the consistent RWA treatment of the 
connecting leads in the presence of ac fields. 

The double-dot model. — We consider the setup 
sketched in Fig. 1 described by the time-dependent 
Hamiltonian H{t) = -ffdots(0 + ^^leads + -^contacts, where 
the different contributions correspond to the quantum 
dots, the leads, and the tunneling coupling to the respec- 
tive lead. We disregard interaction and spin effects and 
assume that intra-dot excitations do not play a role such 
that each dot is well described by a single energy level. 



Then, the double quantum dot Hamiltonian reads 

i^dotsW = (clc2 + 4ci) + ^(clci - 4C2), (1) 

where the fermion operators ci 2 and c\ 2 annihilate and 
create an electron in the left and the right dot, respec- 
tively. The on-site energy difference e(<) = eo + Acos{ilt) 
is determined by the static internal bias eo, the driving 
amplitude A, and the frequency Q. Typical driving fre- 
quencies range up to 100 GHz [11] such that the wave- 
length exceeds the size of the setup and, thus, the im- 
plicitly assumed dipole approximation is well justified. 

The leads are modeled as ideal electron gases, TJicads — 
J2q ^^(cLqCL? + CpqCRq), whcrc cj^ crcates an electron in 
lead ^ = L, R. The tunneling Hamiltonian 
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Cl + Vr„ 



H.c. 



(2) 



establishes the contact between the dot levels and the 
respective lead. Below, we shall assume within a so- 
termed wide-band limit that the coupling strengths F^ = 
2TrJ2q |V«q|^(5(e — eq), £ = L,R, are energy independent. 
To specify the dynamics, we choose as an initial condition 
for the lead electrons a grand canonical ensemble at tem- 
perature T and chemical potentials ^l.r- The influence 
of lead £ is fully determined by the lesser Green functions 



fJ-R = 



FIG. 1: Level structure of the asymmetric double quantum 
dot in a pump configuration. The solid lines mark the relevant 
levels |1) and |2) with the energies ±eo/2. The arrows indicate 
the dominating scattering process. 
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g£g{t,t') = {i /K){c\^(t')ciq(t)) and the tunnel matrix ele- 
ments Viq [15]. More precisely, these quantities enter the 
expressions for the current and the noise in form of the 
correlation function 

- r)m) = dee-'^^'^m (3) 

of the fermionic noise operator = X^g ^gC^9> 

where /^(e) = (1 + exp[(e — ^i()/kBT\)~^ denotes the 
Fermi function [16]. Henceforth, we consider the case of 
zero bias voltage with both chemical potentials located 
midway between the dot levels ±eo/2, i.e., /^l = A*r = 0. 

Resonant electron pumping. — For harmonic driving, 
the Hamiltonian H{t) obeys time-reversal symmetry and, 
hence, each individual scattering process has a time- 
reversed partner which occurs with the same probability. 
Thus, it is tempting to conclude that the net current of 
both partners and, consequently, the pump current van- 
ishes. This, however, is not the case because the driving 
enables energy non-conserving scattering. In particular, 
there exist processes like the one sketched in Fig. 1: With 
the leads initially at equilibrium, an electron from the 
right lead with energy below the Fermi surface is scat- 
tered into a state in the left lead with energy above the 
Fermi surface. This process contributes to the current. 
By contrast, the time-reversed process does not transport 
an electron because the respective initial state is not oc- 
cupied. The net effect is transport of electrons from the 
lower level to the higher level, i.e., from right to left. 

None the less, the pump might vanish due to the pres- 
ence of an additional symmetry, such as generalized par- 
ity {x,t) {—Xjt + n/il.) which relates two scattering 
processes with identical initial energies. Their contribu- 
tions to the current cancel each other [16]. With equally 
strong coupling to the leads, Fl = Fr = F, generalized 
parity is satisfied for H{t) at zero internal bias eo = 0. 
For finite bias eq 7^ 0, however, this symmetry is broken 
and, consequently a finite pump current emerges. More- 
over, this pump current exhibits resonance peaks includ- 
ing higher-order resonances [17]. 

Within our analytical approach, we focus on strongly 
biased situations, eo » A, and driving frequencies close 
to the internal resonances of the double dot, nhfl = 
(cq + A^)^/^ « eo- In this regime, the dynamics of the 
dot electrons is dominated by the second term of the 
Hamiltonian (1) while the tunneling contribution, which 
is proportional to A, represents a perturbation. Con- 
sequently, a proper interaction picture is defined by the 
transformation U{t) = exp[— |(c|ci — C2C2)(/>(i)] with the 
time-dependent phase 



(j){t) = nQt + — sm{Qt). 



(4) 



This yields the double-dot interaction-picture Hamilto- 
nian i/dots(t) = UHt)Ha,^sit)U{t) - inW{t)U{i). The 
transformation U(t) has been constructed such that H{t) 



obeys the time-periodicity of the original Hamiltonian (1) 
while all its other energy scales are significantly smaller 
than hVt. Thus, we can separate time scales and replace 
Hdots{i) within a rotating- wave approximation (RWA) by 
its time average 

-^dots = -^^(ciC2 + 4ci) - - 4^2) (5) 

with 5 = nUQ—eo and the effective tunnel matrix element 

Aeff = {-irJn{A/nn)A, (6) 

where J„ is the nth order Bessel function of the first kind. 

While the lead Hamiltonian is unaffected by the trans- 
formation U{t), the tunneling Hamiltonian acquires a 
time-dependence, .ffcontacts(t) = 5], ^L94gCie~''^(*^/^ 

^R'j'^Rq'^ze*"^^*-'/^. The lead elimination along the Hues of 
Ref. [15], but here for a time-dependent contact Hamil- 
tonian, reveals that the influence of the leads is no 
longer determined by the noise operators but rather 
by Vh/R(t) = e±*'^(*)/2^L/R(i)- Its correlation fimc- 
tion (?7Lyj^(i — T)??L/R(i)) depends not only on the time- 
difference T, but also explicitly on t. The latter time- 
dependence is 27r/r2-periodic and is therefore much faster 
than all other time scales. Hence, we can replace within 
the RWA the correlation function of ?7l,r by its t-average 



(7) 



where 



/L/R,eff(e) 
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can be interpreted as an effective electron occupation 
number of the levels in lead £. At zero temperature, it 
exhibits steps at e = + {k ^ n/2)Ufl and is constant 
elsewhere. 

The RWA provides a mapping of the originally time- 
dependent transport problem to a static one with renor- 
malized parameters. This problem, in turn, can be solved 
by standard procedures: Both the current and the noise 
power can be expressed in terms of the transmission prob- 
ability T(e) of an electron with energy e. For a two-level 
system in the wide-band limit, one obtains 



T(e) = r2|Gi2(e)|2 = 



F^A? 



\{2e-iTy 



^eff 



■<52| 



(9) 



Then, the current defined as the change of the charge 
in the, e.g., left lead, is given by the Landauer-like for- 
mula / = {e/2TTK) /der(e)[/L,eff(e) - /R,eff(e)]; the cor- 
responding expression for the noise power reads [18] 



S=^l deT{e)[ /.,eff(e)[l-/.,eff(e)] 



=L,R 



(10) 



+ [1-T(e)][/L,eft(e)-/R,eff(e)]'}. 
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Note that in the presence of a driving field, even at zero 
temperature, the electron occupation ,eff is not a simple 
step function and, thus, also the term in the first line of 
Eq. (10) contributes to the noise power. A convenient 
measure for the relative noise strength is the Fano factor 
F = S/2eI which characterizes the noise with respect to 
the shot noise level given hy S = 2el [18]. 

For the remaining evaluation of the energy integrals, 
it is important to note that the transmission (9) is prac- 
tically zero for > A^g + + S^. Thus, for hn > 
Aeff , F, 6, the effective electron occupation (8) is constant 
in the relevant energy range and can be replaced by its 
value at e = 0. One obtains close to the nth resonance 
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2h (A2g + r2+52)3 



(11) 



(12) 



where A„ - /L,cff(0) - /R,cff(0) = E\k\<n/2 JU^f^^^) 
with |A„| < 1. Quite remarkably, for resonant driving 
{S = 0), the pump current assumes a maximum while 
the noise power S generally assumes a local minimum; 
cf. Fig. 2(a). This results in an even more pronounced 
minimum for the Fano factor. 

Floquet transport theory. — Before developing an op- 
timization strategy, we corroborate our analytical re- 
sults by an exact numerical calculation within Floquet 
transport theory [16]: Starting from the Hcisenberg 
equations of motion for the annihilation operators for 
both the lead and the dot electrons, one eliminates the 
lead operators and thereby obtains for the electrons on 
the dots a reduced set of equations. These arc solved 
with the help of the retarded Green function obeying 
[ihd/dt - n{t) + iT/2]G{t,t') = S{t - t'), where n{t) 
is the single-particle Hamiltonian corresponding to the 
double-dot Hamiltonian (1). The coefficients of the equa- 
tion of motion for G{t,t') are 27r/i7-periodic and, conse- 
quently, its solution can be constructed with the help of 
the Floquet ansatz |V'q(0) = cxp[{—i€a/fi — -fa)t]\(l)a{t))- 
The Floquet states \(pa{t)) obey the eigenvalue equation 
[n{t) - iT/2 - iM/dt]\(f)a{t)) = (ca - ihja)\<Pa{t)). Its 
solution allows to construct the retarded Green func- 

tion G(t,t') = -{^/h)j:,J^^.At)){i^t(t')\eit - t'). Fi- 

nally, one obtains for the pump current a convenient 
Landauer-like expression with an additional sum over 
the sidebands [16, 19]. Since the symmetrized noise 
correlation function S{t,t') = {[I{t), I{t')]+) depends 
explicitly on both times, we characterize the noise by 
the time-average of its zero-frequency component, S = 
{2TT/n) /o "/^^ dt J drSit, t-T). 

Figure 2(a) depicts the numerically evaluated pump 
current and its noise power. For proper cooling, thermal 
excitations do not play a significant role. Therefore, we 
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FIG. 2: (a) Pump current I (dashed line) and its noise power 
S (solid) at fcsT = as a function of the driving frequency for 
coupling strength F = 0.3A, driving amplitude A — 3.7A, and 
internal bias eo = 5A. The dotted lines mark the analytical 
results (11) and (12). Inset: Blow-up of the lower-left corner 
demonstrating that 7 oc O^. (b) Corresponding Fano factor. 



consider zero temperature only. We find that the cur- 
rent exhibits peaks located at the resonance frequencies 

— {eQ + A'^y/^/nh. This agrees well with our analytical 
results (dotted lines), albeit the RWA predicts the loca- 
tion of the current maxima only to zeroth order in A, 
i.e., at the slightly shifted frequencies Q ~ eo/nh. In the 
adiabatic limit, the pump current vanishes proportional 
to fl^. For the chosen parameters, the noise power S pos- 
sesses clear minima, each accompanied by two maxima. 
In the vicinity of the resonance, the noise is considerably 
below the shot noise level 2el; cf. Fig. 2(b). This fea- 
ture is notably pronounced at the first resonance. Far 
from the resonances, the current becomes smaller and 
the Fano factor is close to = 1. The comparison of 
the numerically exact results with the current (11) and 
the noise power (12) [dotted lines in Fig. 2(a)] leads to 
the conclusion, that the RWA predicts both the current 
maxima and the noise minima sufficiently well to employ 
these expressions for a parameter optimization towards 
low-noise pumping. 

Tuning the electron pump. — We have already seen that 
the condition of large current and low noise is met at 
the internal resonances of the biased double-dot setup. 
Thus, we can restrict the search for optimal parameters 
to resonant driving. As a figure of merit for the noise 
strength we employ the Fano factor for (5 = 



pin) _ 



^ J_ _ Xn r^SAlff - F^) 
2eJ(») 2A„ 2 (A2g + r2)2 ' 



(13) 



which is a function of A„ and Acff/F. The second term is 
minimal for Aeff/F = ^5/3, yielding = 1/(2A„) - 



4 



1.0 



0.5 ■ 



J7_ 
32 

0.0 




r = o.iA 



0.0 



0.5 1.0 



1.5 



FIG. 3: Fano factor F at the first resonance for various cou- 
pling strengths F. The exact Floquet calculation (solid lines) 

is compared with the RWA for 5 = (dashed). The inset de- 
picts the minimal Fano factor in dependence of F for eo = 5A 
at the resonance n = v + e^/Ti. The dotted lines mark 
the optimal Fano factor Fopt = 7/32. 



9A„/32. Thus, the optimal Fano factor is assumed for 
A„ = 1 and reads F^^t = 7/32 » 0.219. 

In the following, wc restrict ourselves to the prime 
resonance (n = 1) for which Aeg = J\{A/h^)/\ and 
Ai = J^{A/2%Q) [20]. Then, the value Ai = 1 is as- 
sumed for A = which means Aog = 0; this unfortu- 
nately imphes a vanishing current (11). Therefore, the 
central question is whether it is possible to find a driving 
amplitude providing on the one hand an appreciably large 
pump current, while on the other hand yielding a noise 
level close to Fopt. The numerical results depicted in 
Fig. 3 indeed suggest this possibiUty: The Fano factor is 
close to the optimal value already for a finite ampUtude. 
A closer investigation reveals that the location of the min- 
imum corresponds to Acff = Ji{A/TiVl)/S. = i,/5/3 F, in 
compliance with our analytical considerations. In partic- 
ular, the minimum is shifted towards smaller values of 
A/hVl for weaker coupling F. Moreover, the RWA solu- 
tions (11) and (12) agree very well with the numerically 
exact results, although they slightly underestimate the 
noise. This discrepancy diminishes as (not shown). 

The data also reveal that in the interesting regime, the 
ratio A/UQ. is considerably smaller than 1 and, hence, 
we can employ the approximations Ja{x) « 1 — .t^/4 and 
Ji(a;) « x/2 valid for small arguments. It is now straight- 
forward to obtain to lowest order in A/TiQ. the expressions 
Aeff = A/\/2nn and F^^) = 7/32 + (^A/lQUVlf. For in- 
stance, choosing A = Q.iMl, the noise level lies merely 
5% above Fopt and the condition Aeff = i/5/3r corre- 
sponds to F « O.IA, i.e., to weak dot-lead coupling. This 
estimate is confirmed by the inset in Fig. 3 which, in addi- 
tion, demonstrates that F f« Fopt for F < O.IA. For such 
a small coupling F, interaction-induced electron-electron 
correlations typically play a minor role. 

In the experiment of Ref. [11], a typical inter-dot cou- 



pling is A = 50/ieV. Then, an internal bias eo = 5A 
corresponds to the resonance frequency Q. = 5A/h « 
27r X 60GHz. Tuning the lead coupling to F = O.IA re- 
sults in an optimized pump current of the order 200pA 
with a Fano factor F « 0.23. 

Conclusions. — The analysis of ac-driven, asymmetric 
double quantum dots demonstrates that optimal pump- 
ing is achieved beyond the adiabatic regime. In partic- 
ular, the ideal modus operandi requires a large internal 
bias at resonant driving in combination with a strong 
inter-dot coupling A > lOF. The resulting pump current 
then assumes a maximum while, interestingly enough, 
the (absolute) noise power assumes at the same time a 
minimum such that the Fano factor becomes remarkably 
small. In order to systematically tune the pump into a 
low-noise regime, we have derived analytical expressions 
for both the current and its noise power. The compar- 
ison with the numerically exact solution fully confirms 
the validity of this approach. Our findings convincingly 
suggest that coupled quantum dots are ideal for pumping 
electrons effectively and reliably at a low noise level. 
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